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We calculate the power spectrum and hemispherical asymmetry of isocurvature and curvaton perturbations 
due to a complex field <I> which is evolving along the tachyonic part of its potential. Using a semi-classical 
evolution of initially sub-horizon quantum fluctuations, we compute the power spectrum, mean field and hemi- 
spherical asymmetry as a function of the number of e-foldings of tachyonic growth AN and the tachyonic mass 
term cH 2 . We find that a large hemispherical asymmetry due to the modulation of |S>| can easily be gener- 
£f>^ , ated via the spatial modulation of |<I>| across the horizon, with A|4>|/|<!>| > 0.5 when the observed Universe 

exits the horizon within 10-40 e-foldings of the beginning of tachyonic evolution and c is in the range 0.1-1. 
The spectral index of the isocurvature and curvaton perturbations is generally negative, corresponding to a red 
, power spectrum. Dark matter isocurvature perturbations due to an axion-like curvaton with a large hemispher- 

ical asymmetry may be able to explain the hemispherical asymmetry observed by WMAP and Planck. In this 
case, the red spectrum can additionally suppress the hemispherical asymmetry at small scales, which should 
make it easier to satisfy scale-dependence requirements on the asymmetry from quasar number counts. 
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I. INTRODUCTION 



(N 

Isocurvature perturbations are a natural feature of a range of dark matter and baryogenesis models. A common class of model 
is based on a complex field <t> whose phase determines the density of the isocurvature component, with axion dark matter and 
l— ~~ 1 j Affleck-Dine baryogenesis being two well-known examples. In these models the isocurvature perturbations are associated with 
quantum fluctuations of the phase angle. 

The nature of the isocurvature perturbations depends on the evolution of <t> along its potential and the time at which the 
. observed Universe exits the horizon. In this paper we will explore the case where <t> evolves from |<t>| = along the tachyonic 
part of its potential, where the potential is of the form V(<t>) ~ — cH 2 \<t>\ 2 . In particular, we will compute the mean field in a 
horizon volume, the spectrum of isocurvature perturbations and the mean hemispherical asymmetry across a horizon volume, 
as a function of the time of horizon exit and of c. We will show that a red isocurvature perturbation spectrum and a significant 
I/"") \ hemispherical asymmetry are generic features of this class of isocurvature perturbation model. 

A topical application of these results is to the hemispherical power asymmetry of the CMB, observed initially by WMAP 
|Q]] and recently confirmed by Planck iHIH]. One possible explanation is a dark matter isocurvature perturbation with a large 
hemispherical asymmetry due to the decay of a curvaton [4]. (See also fH-fH]) 1 . We can apply our isocurvature analysis directly 
J> ' to the case of an axion-like curvaton due to a complex field <t>. A modulation across the horizon by A|<t>|/|<I>| > 0.5 is then 
necessary to account for the observed asymmetry (]4J. We will determine the conditions under which such a large modulation 
can be achieved. The red curvaton perturbation spectrum expected in these models will alter the analysis of the model from 
previous studies, which assume a scale-invariant curvaton spectrum [4]. In particular, a red spectrum should allow the model to 
more easily satisfy scale-dependence requirements on the hemispherical asymmetry from quasar number counts [11] and other 
observations 11211 . 

In Section 2 we introduce the model and the method we use for its analysis. In Section 3 we derive the mean <t> field in 
a horizon volume and the power spectrum and hemispherical asymmetry of the isocurvature perturbations. In Section 4 we 
apply our analysis to obtain the hemispherical asymmetry and power spectrum in specific cases. In Section 5 we present our 
conclusions. Some details of the semi-classical method are included in an Appendix. 



'Electronic address: j.mcdonald@lancaster.ac.uk 

1 Recent alternative proposals to account for the asymmetry are discussed in l9l [Toll . 
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II. EVOLUTION OF ISOCURVATURE FLUCTUATIONS IN A TACHYONIC POTENTIAL 



In our model we consider a complex field <t> = -^j<? to be the origin of the isocurvature perturbations. (We will refer to 
<t> as the 'isocurvature field' in the following.) The density of the isocurvature component of the Universe is assumed to be 
proportional to 0" for some power n. For example, n = 2 in the case of axion-like dark matter in the small misalignment angle 
limit, while n = 1 for baryogenesis in the case where the CP-violating phase is proportional to 9, as in the case of Affleck-Dine 
baryogenesis. 

A closely-related example, which is relevant to the WMAP/Planck hemispherical asymmetry, is a dark matter isocurvature 
perturbation due to the decay of a sub-dominant curvaton field [4]. In the case of an axion-like curvaton field, the curvaton is 
proportional to 0. 

The potential of the isocurvature field is 

v(<*>) = - c i/ 2 |4>| 2 +v,, /r (|<i>|) . (i) 

Here c parameterizes the tachyonic mass term. Vnf t contains the terms which determine the minimum of the potential |<5| m /n. 
These terms will not play a role in our analysis. 

We assume that the isocurvature field is initially at |<t>| = 0, for example by c being initially negative. The field then begins 
to evolve from |<t>| = at an initial time which we define to be t = 0, after which c is positive and the field evolves in the 
resulting tachyonic potential. Although we do not single out any particular class of model, we note that this scenario is a natural 
possibility in the context of supergravity. A — cH 2 \<t>\ 2 term is a natural supergravity correction and c can change sign if the 
dominant inflaton field changes at some time, with the later inflaton field having different Kahler couplings to the isocurvature 
field. 

The field O can be expanded in terms of real fields, 4> = ((j)i + ityi)/ \2, Since the evolution is linear during the tachyonic 
regime, it is the same for §\ and (j)2, so we will consider a generic real scalar field (]). The initial state of the field at t = is a 
Bunch-Davies vacuum on subhorizon scales. We will assume that no superhorizon perturbations can form before the tachyonic 
growth era, which is true if |c| > 1 before the onset of tachyonic growth. 

In order to evolve the field from this initial state, we will use a semi-classical method that is commonly used in numerical 
simulations of tachyonic preheating 1U3M16I1 . The method is based on the Wigner function for the field. The Wigner function 
method replaces the quantum field with a classical field which evolves from initial conditions which have a classical probability 
distribution. Quantum correlation functions in the semi-classical limit are then equal to the classical correlation functions 
averaged over the classical initial probability distribution. We review the method in the Appendix. 

The initial classical modes are Gaussian distributed with a random phase. Although the initial conditions are sub-horizon and 
not initially semi-classical, the modes can still be evolved from classical initial conditions, provided that the mode enters the 
semi-classical regime while the field is still undergoing linear evolution i.e. the field is still dominated by the tachyonic mass 
term 11151 Il6ll . This is because the Wigner function in the semi-classical limit becomes a classical phase space density which 
must reproduce the late-time semi-classical behaviour. Therefore the final result for late-time semi-classical modes will be the 
same, regardless of when the classical initial conditions are applied. 

The rms values of the initial classical modes in de Sitter space are (Appendix) 

1 / //2\ 1/2 

l?(M)U« = - 

and 



|?(M)U = ._^ + __] (2) 



\p^)\ rms = — \ [ l + —j , (3) 

where T] = — 1/aH is the conformal time and initially r| = — l/H. 

Since we will be interested in mean quantities, we will evolve the classical modes from the RMS initial conditions Eq. (O and 
Eq. (f3]l, with the initial time defined to be t = and the initial scale-factor defined to be a = 1 . The random phases will not play 
a role in our analysis. 

The classical modes during the tachyonic era are defined by 

(Kx,a) = 47E^(ay k - x , (4) 
VV k 

where we use the scale factor to parameterize time. (We define modes in a box of volume V and side L, taking the continuum 
limit where appropriate.) The mode equation is 

9 2 frk | 49fr k | k 2 c 
da 2 a da a 4 H 2 a 2 
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The general solution is 



(6) 



where J v (x) are Bessel functions. To fix c\ and C2 we match this to the rms initial conditions, Eq. (O and Eq. ||3}, at a = 1. 
Modes corresponding to the observed Universe will be subhorizon at a = 1 . In this case k > 2kH and we can use the asymptotic 
form of the Bessel function, 



J v {x) 



2 ( V7I 7C \ „ / 1 
COS [X - — - - + 0\ -rpr 

Six V 2 4/ Vx 3 / 2 



; x> 1 



Therefore, for k/aH ^> 1, which is generally true for modes at a = 1, 

k V4c + 971 7C 



* k(fl)=cl ^ cos U 



2atf / £ V4c + 9tc 7C 

C2A / COS 

nk \aH 4 4 



(7) 



(8) 



<j)(x,a) = y(x,-l/.ff) and (j)(x,a) = n(x,-l/H) at a = 1 and r| = -1/H. Thus (j) k (l) = |#(k, -l/7/)| rmj and (j>k(l) = 
|p(k, — 1/H)\ rms . Therefore the initial conditions at a = 1 are 



Ml) 



ff 2x 1/2 



(9) 



and 



These determines ci and C2 to be 



and 



where 



3a 



(1) 



C2 = 



AH (c+i_ — s+c_) 

7C (S+-C+¥) 

4H (c+s_ — s+c_) 



. k , V4c + 971 7c\ . / V4c + 97C 7t 

c ± = co S ( -± ^ + -j ; s± = sm (-±^— + - 



and c + i_ — i + c_ = — sin(7C\/4c + 9/2). 



III. ISOCURVATURE PERTURBATION SPECTRUM AND HEMISPHERICAL ASYMMETRY 



A. Mean field 



(10) 



(11) 



(12) 



(13) 



As the modes grow, there will be a mean value of the field in a horizon-sized volume at a given a, due to modes which have 
wavelengths larger than the horizon. The rms field is given by 



<( K x > f ) 2 >= (2^)3 / M 2 d'k 



We define J{k/aH) by 



(14) 



(15) 
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where 



Then 



and 



AH 

C\,2 = \ Ci,2 ■ (16) 

V K 



<(j)(x,f) 2 >= ~ / J 2 [ — )k 2 dk. (18) 

^ j 8%HalJ knUn \aHj 

Here & m ,„ = 2ft// corresponds to initial the horizon-sized mode at a = 1 and k max = 2%aH is the horizon-sized mode at a. 
Changing variable to y = k/aH, we obtain 

f(a)=<^(x,t) 2 >= I ^K l (a), (19) 



where 

c2% 



K l {a) = J u J 2 (y)y 2 dy. 



(20) 



-r2. 



It is useful to have an approximation for (j) (a) which shows the a dependence explicitly. To do this we split the integral into 
two parts, 

K l {a)=K l0 +K l {a), (21) 

where 

[■2% 

Kio= J 2 (y)y 2 dy (22) 

J27ia 

and 

r2na 

Ki(a) = J u J 2 (y)y 2 dy. (23) 



We find that K\ is an essentially a-independent constant for each value of c. a is chosen small enough that the modes in K\ (a) 
satisfy y -C 1 . In this case we can use the small x form for the Bessel functions, 

^> = n^T)©"( 1+0 © 2 )' <24 > 

Therefore 



/(y)«ci , - (j-) +c 2 , - (j-) 2 . (25) 

w r(i - iv 7 ?^) V2^ 2 r(l + iV4?+9) V2^ 

Fory < 1, 

/(y)wci , - (j-) . (26) 

w r(i-±V47+9) w 

Since k/H > 2%, we can approximate ci by 

s- sin (4) 

Cl« ? ; r « ? V g^-N • (27) 

sin(5^f±2) sin(^f±5) 
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Therefore 



Kx{a) 



r 2(l- sin 2 (2^5) 



2jta 



sin 2 (ay)y 2 



dy. 



(28) 



Since sin 2 (ay) is rapidly varying when a 3> 1, we can replace it with its average value sin 2 (ay) = 1 /2. With this we find 



*i(a) 



/4c+9-l 



/4c+9-3 



.2jtJ 



r 2 (l-iV4?+9)sin 2 (2^) (\/4c + 9-3) 
Thus an approximation showing the a dependence explicitly is 



(29) 



Ki(a) = A+Ba y 



■1) 



71 



(30) 



where A = ^1 , y = y/4c + 9 — 3 and B is the a-independent factor in Eq. (l29l . We give the values of A, B, and y as a function 
of c in Table 1. Using a = 0.02, we find excellent agreement with the exact value of K\ (a) for a > 50. (In practice we will be 
interested in a 3> 50.) 



c 


A 


B 


Y 


0.1 


7.64 


4.49 


0.066 


0.2 


7.77 


2.12 


0.130 


0.3 


7.94 


1.34 


0.194 


0.5 


8.26 


0.72 


0.317 


0.75 


8.78 


0.45 


0.464 


1.0 


9.51 


0.29 


0.606 



TABLE I: Parameters for the approximation to K\ (a) when a = 0.02. 



B. Power Spectrum 

At a given value of a, the (j)i and (j)2 fields will evolve to mean values of the order of (j) in a given horizon volume. We can 
then make a field redefinition so that (j)i has non-zero mean value (the radial direction) while (j)2 has only fluctuations about zero 
mean (the angular direction). Then 



<j)i «<j) + 8(|)i 



(31) 



and 



(|>2 » §<j)2 



(32) 



The angular field fluctuation is then 



802 80 



59 = w ^ h-ni) 



+ 801 



80i (x) 







(33) 



Thus 88 = 802/0 will give the mean isocurvature perturbation while the spatial modulation of the isocurvature perturbation is 
due to 80i (x), 



A(89)(x) 




80! (X) 


86 








(34) 
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The isocurvature perturbations are generated when 8(j)2 fluctuations exit the horizon. Thus the power spectrum is determined 
by the superhorizon fluctuations of 8(j)2. For superhorizon modes we can use the small x expansion of the Bessel function as 
before. Therefore 



< 8$ > = 



1 



16jti/fl 3 r 2 (l - !V4F+9)sin 2 



2aH 



k 2 dk< 



dk 



(35) 



Therefore the spectral index of the isocurvature perturbation is 



ritso = 4 - V47+9 . 



(36) 



This is less than 1 for all c > 0, corresponding to a red isocurvature power spectrum. We give the spectral index as a function of 
c in Table 2. 



c 


0.1 


0.3 


0.5 


0.75 


1.0 


n iso 


0.934 


0.806 


0.683 


0.536 


0.394 



TABLE II: Spectral index of the isocurvature perturbation for c in the range 0.1 to 1.0. 



C. Hemispherical Asymmetry 

The hemispherical asymmetry of the isocurvature perturbation is due to the difference in 88 across the horizon at present. 
This is due to the spatial change in 8(j)i /(j) across the horizon when the observed Universe exits the horizon. To make this clear, 
consider the evolution of 8(j)i /(j) once the 8(j)i fluctuation has wavelength larger than the horizon. Within a horizon volume, 8(j)i /(j) 
is approximately spatially constant and, since the evolution is linear during the tachyonic regime, this ratio remains constant as 
(j) grows. Therefore the modulation over the volume of the present Universe is the same for all subsequent 88 perturbations as 
they exit the horizon and become fixed in amplitude. 

We therefore calculate the rms value of 8(j)i(x + 8x) — 8(j) i (x) for |8x| corresponding to the comoving horizon (aH)^ 1 at a 
given value of a. The mean squared value is 

< (8(j)i(x + 8x) -8(j)i(x)) 2 >=< 8(j)i(x + 8x) 2 > + < 8(j)i(x) 2 > -2 < 8(j)i(x + 8x)8(^i(x) > 



= 2[<8(j)i(x) 2 >-<8(j)i(x + 8x)8(j)i(x)>] . (37) 

< 8(j) 2 (x) > is given by 

< Stf(x) >= ^ I \^ k \ 2 d 3 k . (38) 

The cross-correlation is given by 

< 8^) 1 (x + 8x)8^ 1 (x) >= ^ j \$ k \ 2 e ikM d 3 k . (39) 

Therefore 

<(8^ 1 (x + 8x)-8^) 1 (x)) 2 >-^3 J \$ k \ 2 (l- e ' kM )d 3 k=-^ J H) k | 2 (l-cos(k.8x))A. (40) 

Sx is along the radius of a horizon volume centred at 8x = 0, with |8x| = (aH)^ 1 . For different k directions with equal value of 
k, we can average over the directions, therefore we make the mild approximation that k.8x = A:|8x| cos(o) « A:|8x|/2, where o is 
the angle between k and 8x. Therefore, integrating over superhorizon modes, we obtain 

< (54>i(x + 5x) -8^(x)) 2 >= ^ W 2 sin 2 f-M k 2 dk . (41) 

JKnin \2aHJ 
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Using Eq. ( TTTb and changing variable to y — k/aH, we find 

Sjf 2 =< (8^) 1 (x + 8x)-8^ 1 (x)) 2 >= ^K 2 (a) , (42) 

where 

£ 2 (a) = J^J 2 (y) S m 2 {^)y 2 dy . (43) 

We find that /^(a) is constant with respect to a for values of a > 10. Thus the spatial modulation of the field across a horizon 
radius becomes independent of time, whereas the mean field increases with time. In Table 4 we give the values of K 2 for c in the 
range 0.1 to 1 . 



c 


0.1 


0.2 


0.3 


0.5 


0.75 


1.0 


K 2 


3.56 


3.58 


3.59 


3.65 


3.72 


3.80 



TABLE III : K 2 for values of c in the range 0. 1 to 1 .0. 



IV. RESULTS 

There are two cases of interest: (i) the possibility of observing isocurvature perturbations with a red spectrum and a significant 
hemispherical asymmetry 2 and (ii) the specific case where a dark matter isocurvature perturbation with a large hemispherical 
asymmetry due to curvaton decay may account for the hemispherical power asymmetry observed by WMAP and Planck. 

A. Isocurvature spectrum and hemispherical asymmetry 

In Figure 1 we show the rms spatial modulation of (j) across the horizon, A(j)/(j), as a function of the number of e-foldings AN 
after the beginning of tachyonic evolution, for c in the range 0. 1 to 1 . In Tables 4 and 5 we give the value of AN as a function of 
c below which a modulation A(j)/(j) equal to 0. 1 and 0.5 can be achieved. 

Rather generally, we find that a significant hemispherical asymmetry of the isocurvature perturbation can be achieved. Values 
of A(|)/(j) > 0.1, corresponding to a greater than 10% hemispherical isocurvature perturbation asymmetry across the horizon 
radius, are obtained when AN < 100 for c = 0.1 and AN < 15 for c — 1.0. Therefore if our Universe exits the horizon when AN 
is within these limits, a greater than 10 % asymmetry in the isocurvature perturbation can be achieved. This asymmetry would be 
associated with a red isocurvature spectrum, therefore the combination of a hemispherical asymmetry and red spectrum would 
be a signature of this class of model. 

In Figure 2 we show the value of (j)/// as a function of AN for c in the range 0.1 to 1. In general, the value of (j)/// at which 
large A(j)/(j) occurs is not larger than 7. Specifically, for A(j)/(|) = 0.1 we find that (j)/// is approximately 7 for all c. Thus even 
though a large number of e-foldings of tachyonic growth can occur, the change in the mean field is relatively small. Therefore is 
quite plausible that the isocurvature field will not reach the minimum of its potential before the large asymmetry is generated. 

B. WMAP/Planck hemispherical asymmetry 

In J3l it was proposed that the WMAP hemispherical asymmetry might be explained by a dark matter isocurvature perturbation 
due to a curvaton field with a large hemispherical asymmetry. In the case of an axion-like curvaton, this would require a 
modulation of the curvaton across the horizon radius by A(|)/(j) > 0.5. From Figure 1 and Table 5, we see that it is possible to 
generate large modulations of the curvaton in this model. For c = 0.1, A(|)/(j) > 0.5 is generated if AN < 40, while for c = 1 this 
is achieved if AN < 10. 



2 We calculate the hemispherical asymmetry in the rms magnitude of the isocurvature perturbation. The hemispherical isocurvature power asymmetry will then 
be a factor of 2 larger, since the power is proportional to the square of the perturbation. 
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A strong constraint is from quasar number counts, which show that the hemispherical asymmetry must decrease at large / 
fiUl . This can be achieved if the asymmetry is due to an isocurvature perturbation, since the isocurvature power naturally decays 
relative to the adiabatic power at large / fl. In the curvaton model for dark matter isocurvature perturbations, the curvaton 
decays after dark matter freeze-out, producing equal adiabatic and isocurvature contributions to the total power spectrum, both 
of which have a hemispherical asymmetry. This leads to a strong constraint on the possible hemispherical asymmetry from the 
quasar bound, since the adiabatic asymmetry does not decrease with scale and must satisfy the quasar bound at large /. 

However, the analysis of [4] assumes a scale-invariant curvaton spectrum. The effect of a red spectrum will be to decrease 
both the isocurvature and adiabatic power on larger scales. This should therefore relax the scale-dependence bound on the total 
hemispherical asymmetry, allowing a larger hemispherical asymmetry and so a smaller contribution to the total CMB power by 
the dark matter isocurvature perturbation. 




20 40 60 80 100 120 140 160 

A N 

FIG. 1: Values of A§/§ versus AN for c is the range 0.1 to 1.0. 



c 


AN 




0.1 


83.1 


6.84 


0.3 


34.6 


6.86 


0.5 


23.2 


6.94 


0.75 


17.1 


7.06 


1.0 


13.8 


7.10 



TABLE IV: Number of e-foldings of tachyonic growth and the mean field when A<])/<]) = 0.1. 



V. CONCLUSIONS 

We have studied the properties of the isocurvature perturbation generated by a complex field which grows from <t> = in a 
tachyonic potential. The model has one free parameter, the tachyonic mass term, cH 2 . 

A general feature of this class of model is that the spectral index of the isocurvature perturbations is less than one, with the 
spectrum becoming increasingly red as c increases, with «,■„, = 0.934 for c — - 0.1 and «,■„, = 0.394 for c = 1 .0. 

We find that a hemispherical isocurvature perturbation asymmetry is a natural feature of isocurvature perturbations in this 
model. Asymmetries of the isocurvature perturbation greater than 10% across the horizon can be easily achieved if the observed 
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FIG. 2: Values of §/H versus AN for c in the range 0.1 to 1.0. 



c 


AN 


iSf/H 


0.1 


36.8 


1.50 


0.3 


18.6 


1.52 


0.5 


13.3 


1.52 


0.75 


10.3 


1.54 


1.0 


8.6 


1.54 



TABLE V: Number of e-foldings of tachyonic growth and the mean field when A<|>/(]) = 0.5. 



Universe exits the horizon during inflation at between 14 and 80 e-foldings after the beginning of tachyonic evolution for c in 
the range 0. 1 to 1 .0. 

Therefore a signature of this class of isocurvature perturbation model would be the observation of a red isocurvature spectrum 
together with a hemispherical isocurvature perturbation asymmetry. 

A topical application of this model is to the hemispherical power asymmetry of the CMB, first observed by WMAP and 
recently confirmed by Planck. It has been proposed that this might be explained by a dark matter isocurvature perturbation due 
to the decay of a curvaton which has a spatial modulation across the present horizon radius of at least 50% fl. 

Our analysis of isocurvature perturbations can be directly applied to the case of an axion-like curvaton field. We find that 
large hemispherical asymmetries of the dark matter isocurvature perturbation can easily be achieved. For the case c = 0.1, a 
50% modulation of the curvaton field requires that the observed Universe exits the horizon no later than 37 e-foldings after the 
beginning of tachyonic evolution, while for c = 1 .0 it must exit before 9 e-foldings. 

A possible advantage of the tachyonic model is that the red spectrum can additionally suppress the isocurvature power and 
hemispherical asymmetry at large /. In the analysis of [4], an important constraint comes from quasar number counts and other 
observations, which show that the hemispherical power asymmetry is smaller at larger /. In the case where the dark matter 
isocurvature perturbation comes from curvaton decay, there is an equal adiabatic perturbation which also has a hemispherical 
asymmetry. Since, unlike the isocurvature perturbation, this does not decay at large I, the quasar bound limits the amount 
of adiabatic hemispherical asymmetry (and so dark matter isocurvature asymmetry) for a given adiabatic perturbation from 
curvaton decay [4]. The additional suppression due to the red spectrum should therefore relax the quasar bound and allow a 
larger hemispherical asymmetry. This could in turn relax constraints from the upper bound on the isocurvature perturbation. 

It remains to be established whether an isocurvature perturbation from curvaton decay can consistently account for the 
WMAP/Planck hemispherical asymmetry while satisfying the latest Planck constraints on the isocurvature perturbation JH. 
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Our model shows that a large hemispherical asymmetry in the isocurvature perturbation can be generated in a plausible model. It 
also illustrates the need for any future analysis to include a possible red spectrum for the isocurvature and curvaton perturbations. 
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Appendix. Wigner function method for semi-classical tachyonic mode evolution 

We use the same method used in numerical studies of tachyonic preheating and oscillon formatio n lfl3l [T5[ [l6ll . The original 
formulation is given in 01711. A summary of the methods is contained in itlal and in the Appendix of IU5I1 . 

For consistency with lll5l UM . we will consider a finite box of volume V and side L, taking the continuum limit where 
appropriate. The modes are defined using periodic boundary conditions, with wavenumbers ki = 2%rii/L. The field <t> can be 
expanded in terms of real fields, <I> = (<|)i + ityi)/^/2. Since the evolution is linear during the tachyonic regime, it is the same for 
(|)i and (|)2, so we will consider a generic real scalar field (]). The action is 



S = 



d 4 X\ 



(A-l) 



This can be expressed in terms of y — a§ and integrated over comoving coordinates and conformal time 



d x 



lli a 



cH z a 



2„2 



(A-2) 



where primes denote derivatives with respect to T|. 

This is canonically quantized in terms of y and its conjugate momentum K = dL/dy. The mode expansion of the fields is 



y(k,T 1 ) = ^=^ ? (k, 11 y k - x 
vv k 



(A-3) 



and 



W k 



(A-4) 



where 



pOMi) = « (M)--9(M) 



(A-5) 



These are quantized via 



«(M) =/kak+/ka[ k 



(A-6) 



and 



P(k^) = -igkak + igta'_ k 



(A-7) 



where [ak,fl] [] = §k,i- The mode functions satisfy the classical field equation 



ji a 
a 



/k--/k = -(k 2 + mV)/ k 



(A-8) 
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In the case of de Sitter space with constant H, and defining T| = — 1 /aH, this becomes 



A 



2-- 



A = o 



(A-9) 



The correctly normalized solution for quantum modes in de Sitter space is then 

A(n) = ^e i (f v +?)y=nff«(-)bi) ) 



(A-10) 



where v 2 = 9/4 — m 2 /H 2 . In the following we will consider the field to be effectively massless on subhorizon scales prior to the 
tachyonic growth era, therefore 



£ 2 y/2k 



1 



(A-ll) 



and 



-ikr\ 



(A- 12) 



The Wigner function method replaces the quantum modes by classical modes at an initial time. The Wigner function is a 
quasi-probability distribution for quantum states, W(q,p,t), which gives the probability for q on integrating over p and vice- 
versa. In the semi-classical limit the Wigner function W(q,p,t) becomes a classical phase space distribution. However, we 
can define the classical phase distribution for the late-time semi-classical state at any time during the linear evolution of the 
fields, even when the semi-classical limit is not valid. This is because taking/; — > does not affect the late-time semi-classical 
state. Therefore if we take h — > at earlier times, the resulting classical phase space distribution must reproduce the correct 
semi-classical limit at late times. The only requirement is that the quantum fields must enter the semi-classical regime during 
the linear evolution of the fields, since it is only during linear evolution that we can evaluate the classical probability distribution 
from the Wigner function. We also note that this is an exact solution of the field theory in the semi-classical limit; there is no use 
of perturbation theory. 

The resulting classical modes are Gaussian distributed with a random phase. The distribution is given by MM 



1 

P{qk,t) = -j exp 



Ml 
iaP 



(A- 13) 



with pk = fk^k/|/k| 2 and — Im{f^gk) — 1 /2kx\. Thus the rms classical modes at conformal time r| are 



k(M)U 



1 v 1 



V2k V k W 



1/2 



(A- 14) 



and 



b(k,r|)|™ 



1 f l- 1 



2£|r|| V k W 



-1/2 



(A-15) 



These give the initial conditions at t = 0, corresponding to r| = — l/H, for the classical evolution of the modes in the tachyonic 
era. 
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